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Abstract. 

Neutron structure functions can be extracted from proton and deuteron data 
and a representation of the deuteron structure. This procedure does not require 
DIS approximations or quark structure assumptions. We find that the results 
depend critically on properly accounting for the Q 2 dependence of proton and 
deuteron data. We interpolate the data to fixed Q 2 , and extract the ratio 
of neutron to proton structure functions. The extracted ratio decreases with 
increasing x, up to x m 0.9, while there are no data available to constrain the 
behavior at larger x. 
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1. Introduction 

Traditionally, much of our information regarding the quark structure of the nucleon 
has been gleaned from charged lepton scattering. If the one-photon exchange 
approximation is sufficient, then inclusive structure functions of a target nucleon 
or nucleus are well defined invariant functions of the cross sections. The structure 
functions are dimcnsionless invariant functions of the target four-momentum P and 
the four-momentum transfer q. Thus nucleon structure functions are functions of 
two invariants, e.g. the momentum transfer, Q 2 = —q 2 , and the energy transfer, 
v = P ■ q/m, where m = VP 2 is the nucleon mass. The structure functions can also 
be taken as functions of other invariants, e.g. Q 2 and x, where x = Q 2 /2P ■ q. It 
follows from the definition of x that 

independent of any structure assumptions. For "deep-inelastic" scattering, where 

£ = £^ = -£_»l (2) 
Q 2 ~ m 2 Q 2 ~ Am 2 x 2 ' V ' 

the parton model implies that the structure functions depend only on x ( "scaling" ) , 
with a weak logarithmic dependence on Q 2 ( "scaling violations" ) due to QCD evolution 
of the parton distributions. 

Since the discovery of scaling in deep inelastic electron-proton scattering [IJ [2j [3l 
IHE], there have been a number of measurements 6, 7, 8, 9, 10 of the proton structure 
function. The relation of the proton structure function to the quark structure of the 
nucleon is model dependent. At medium energies and large x, e.g. Q 2 = 10 GeV 2 and 
x = 0.9, yielding Q 2 /4m 2 x 2 = 3.5, the limit of © is not satisfied and the scattering 
is not "deep inelastic" . 

To test models of the quark structure of nucleons it is important to have 
experimental results for both protons and neutrons. The nucleon structure functions 
for x > 0.4 are particularly important for constraining models. Since neutron targets 
are not practical, experimental neutron structure functions must be extracted from 
deuteron and proton scattering data, without any quark assumptions. This requires 
consistent Poincare covariant representations of both the deuteron states and the 
current-density operators. 

Realistic nucleon-nucleon potentials are derived from Lagrangians which 
determine two-nucleon Green functions, and the interaction dependent exchange 
currents [111 112) . Different truncations, regularizations and Lorentz invariant 
constraints are used for the construction of effective potentials. Parameters are 
adjusted to accurately fit both deuteron properties and nucleon-nucleon scattering 
data. These potentials define two-nucleon rest-energy operators (mass operators), 
which are both Lorentz invariant and Galilean invariant. The difference between 
the binding energy Mj_ — 2m and (Mj — 4m 2 ) /4m is a negligible relativistic effect. 
Spin operators are invariant under Galilean boosts, and undergo Wigner rotations 
under Lorentz boosts. This difference is the source of important relativistic effects. It 
does not, however, affect inclusive structure functions of unpolarized targets. If these 
mass operators fit the same scattering and bound state data then they are related 
by unitary transformations. Such unitary transformations modify the separation of 
current density operators into one- and two-body operators. 
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Given the dynamics specified by an invariant mass operator, unitary 
representations of the full Poincare group obtain readily by specification of the 
generators as functions of the mass operator. The generators of a kinematic subgroup 
are independent of the mass operator. This choice of kinematics which affects 
the representation of the current density operators, can be exploited to simplify 
the relations of deuteron structure functions to nucleon structure functions. The 
identification of single- nucleon currents is representation dependent. With null- plane 
kinematics [13] the representations of Lorentz transformations that leave a null vector 
n = {1, — ft} invariant are independent of the dynamics and only the momentum 
component P~ = I ■ P, I — {l,ft} depends on the mass operator. The null vector 
n can be chosen such that Q + = n ■ q = 0. Thus Q 2 = Q\ — Q + Q~ = Q\ is 
independent of the dynamics. Kinematic Lorentz transformations which change P± 
do not affect Q + = 0. With this form of kinematics the impulse assumption for the 
relevant components of the current density is consistent with current conservation and 
Lorentz covariance. 

This form of kinematics permits a Lorentz invariant convolution relation of the 
nucleon and deuteron structure functions, which does not require the deep-inelastic 
approximation, and is independent of the mass spectrum of the final state. Such 
"smearing" is Q 2 dependent but converges to the familiar Q 2 independent convolution 
when Q 2 jv 2 is negligible. It is convenient to choose P + = 2m and P± such that 

P-q=-P±-Q±, v = -^-Q±. (3) 

TO 

The values of invariants are, of course, independent of such a choice. 

In this paper we extract the neutron structure function from proton and deuteron 
data taking the impulse assumption and an effective two-nucleon mass operator. The 
main focus is on the consequences of systematic uncertainties and the Q 2 dependence in 
the proton and deuteron data. A systematic exploration of the uncertainties associated 
with using different models of deuteron structure is beyond the scope of this paper. 

2. Cross Sections and Inclusive Structure Functions 

For inclusive electron scattering from an unpolarized target nucleus the Lorentz 
invariant inclusive cross section can, in the one-photon exchange approximation, be 
expressed in the form 

_^ = 4a^Wfl/2 

dQdE' Q A M Wtfuu, W 

where M — y/P 2 is the target mass and the electron current tensor, 
_ 1 f kfikl + k'^k v \ 
7m* - 2 ^ {k . fc ,j , W 

is a bilinear function of the initial and final electron momenta k and k' . It follows 
from Lorentz covariance and current conservation that the current tensor F^ (P, q) is 
a function of two invariant functions, F\{v, Q 2 ) and F^iy, Q 2 ), 

/n^n v \ P^ P v 

F^(P, q) = - g^)F x {y, Q 2 ) + — F 2 (u, Q 2 ) , (6) 

V q z / P ■ q 



where 



P ■ a 

P = P-^q. (7) 
q 2 
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In the frame specified by ^ and Q± = {\/Q^,0}, the invariant functions are 
proportional to single components F 22 and F ++ = n tl n v F ,u ' of the current tensor 

F++(P,q) 



F 1 (v,Q*)=F M {P,q), F 2 (v,Q A ) = (P ■ q)- 



(8) 



Since q^ /„„ = it follows from these definitions that 



da Aa 2 E' 2 M cos 2 9/2 



dQ.dE 1 



Q 4 (P ■ q) 



F 2 {^Q 2 ) 



2P-q tan 2 6/2 
M 2 



Fi(v,Q 2 ) 



and with the ratio of longitudinal and transverse cross sections R = ol/ot, 



1 + R(u, Q 2 ) = 



M 2 F 2 {u,Q 2 ) 



da 



P-qF x {v, Q 2 ) 
4a 2 E' 2 M cos 2 6/2 



dttdE' 



Q+P-q 



Fo 



i + 2 } Q2+ ; 2) 2 t^6/2 

(1 + R)Q 2 1 . 



(9) 



(10) 



(11) 



The relation of structure functions to cross sections does not involve any assumptions 
about the target structure. 



3. Proton and Deuteron Data 



To extract neutron structure functions using convolution relations we need proton 
and deuteron structure functions at a fixed value of Q 2 . To minimize systematic 
uncertainties, we use a fit to the proton structure function, F 2p , from a recent 
global analysis and direct measurements of the ratios of deuteron to proton structure 
function, Rd p = F 2 d/F 2p . Therefore, only the ratios need to be interpolated to 
fixed values of Q 2 , minimizing the size of these correction. In addition, we minimize 
sensitivity to normalization uncertainties between data sets by using only measurement 
of F 2c i and F 2p taken from the same experiment when forming the deuteron-to-proton 
ratios. 

For F 2p (x,Q 2 ), we use the fit from [2], evaluated at Ql = 12 GeV 2 , which 
provides parameterizations for F lp and F 2p . The parameters are fit to a large body 
of data P3 E OH OH El HI]. The Q 2 dependence comes from the fit to the data, 
rather than from a specific model of scaling violations. We compared the fit to the 
F 2p extraction from several SLAC measurements [El [19] and BCDMS [8] . For both 
experiments, the data are taken from the analysis of Whitlow [H)l [T5] . where a relative 
normalization factor between different experiments is determined and all data are 
extracted using a common parameterization for R = a^/ax- The fit is consistent with 
the data within the quoted uncertainties, yielding xt — 0.86 when compared to all 
data for 6 < Q 2 < 40 GeV 2 . 

For the ratio Rd P (x, Q 2 ), we take the measurements from SLAC [TB], BCDMS [5J 
[20] . and NMC [?]■ The SLAC data are from a reanalysis of several experiments, 
and included a fit to the relative normalizations for the different SLAC experiments. 
While the SLAC, BCDMS, and NMC used different models for R = a L /a T , this has 
no impact on the structure function ratios, as all of the extractions used identical 
values for R is the proton and deuteron, as suggested by world's data [TUl [15] . 

In combining all of the Rd P (x,Q 2 ) measurements, we have to account for 
both statistical uncertainties and the relative normalization of the data sets. The 
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uncertainties are broken up into statistical, correlated systematic, and normalization 
uncertainties. We take the statistical uncertainties to be the combination of the 
counting statistics and the uncorrelated systematic uncertainties. The correlated 
systematic uncertainties come mainly from the relative normalization of the different 
data sets, which yields uncertainties that are typically highly correlated between 
neighboring x points. The detailed analysis [19] yielded a 1% overall normalization 
unceratinty, with smaller relative normalization uncertainties. In combining the results 
with the BCDMS and SLAC measurements, we also fit a relative normalization offset 
between the experiments. In the final analysis, we scale the SLAC ratios down by 0.8%, 
and the BCDMS and NMC results up by 1.9% and 0.1% respectively, in all cases within 
the systematic uncertainties quoted for the ratios. We take a 1% overall normalization 
uncertainty on the ratio, and use uncertainties in the relative normalizations of 0.5% 
for NMC, 1.0% for BCDMS, and 0.4-1.1% for the SLAC experiments, as determined 
in Ref. [H]. 

Because the proton structure function parameterization does not include a 
contribution associated with elastic scattering, the contribution from quasielastic 
scattering in the deuteron would not be reproduced. The kinematics are chosen such 
that this contribution is extremely small. A fit to quasielastic scattering data, made 
to estimate quasielastic and inelastic contributions in the analysis of large x inclusive 
scattering data [2T], indicates that the quasielastic contribution to the Rd P ratio is 
negligible, with a contribution that is at most 0.3% (for the largest x values), and is 
always much less than the experimental uncertainties. 

There are several points in the results of [19] where the quoted values for F 2 d and 
F 2p do not match the value given for the ratio F 2 d/F 2p . This usually corresponds to 
cases where the x or Q 2 values of the hydrogen and deuterium data did not match 
exactly, and corrections had to be applied to form ratios at identical kinematics. This 
correction can be large, especially at large x where F 2p drops rapidly with increasing 
x. This correction is model-dependent, and can introduce a significant uncertainty if 
the correction is large. In addition, there are a few cases where the discrepancy is too 
large to be explained by the interpolation (up to 25% in one case) , suggesting an error 
in either the quoted structure functions or the ratio. To avoid introducing possible 
bias due to such errors and to minimize the model-dependent corrections associated 
with large interpolations, we eliminated points where the quoted value for the ratio 
is more than 2% from the value calculated from the quoted F 2 d and F 2p values. The 
impact of this cut is negligible except at very large x, where the effect is still small 
compared to the experimental uncertainties. As the final result is insensitive to the 
exact value of the cut, we do not include any additional uncertainty associated with 
this cut. 

This procedure gives us a large number of individual values of the ratio Rd P (x, Q 2 ), 
which need to be interpolated to a fixed Q 2 value, Q 2 ,, and binned in x. We choose 
Qq = 12 GeV 2 , as this is the average Q 2 value of the large- ir data, and limit the overall 
Q 2 range so that the average Q 2 in a given bin is within a factor of two of Qq, to 
minimize the size of the interpolation correction. 

It has been observed that the Q 2 dependences of the proton and nuclear structure 
functions [22l [23] EU [25] and structure function ratios [26] are much smaller when 
taking fixed values of £ = 2x/(l + yl + Am 2 x 2 /Q 2 ) rather than fixed x. Thus, 
we fit the Q 2 dependence at fixed £ over the full Q 2 range of the data, from 3- 
230 GeV 2 , to extract the Q 2 dependence. We find dR dp / d{\n{Q 2 )) w -0.013, with 
little x dependence, except for x & 0.1 and x ^ 0.7, where the Q 2 dependence appears 
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to be smaller, but is not well measured. This Q 2 dependence is used to interpolate 
the individual data points to Q 2 y 

For the final analysis, we use only data with W 2 > 3 GeV 2 and 6 < Q 2 < 40 GeV 2 , 
yielding a maximum interpolation correction of <1%, and an RMS correction of 0.4%. 
In the regions where the Q 2 -dependence is not well measured, this correction is small 
compared to the uncertainties, and the associated uncertainty is negligible. After all 
of the cuts, we have 450 Rd P measurements, all corrected to Q 2 ,. We then bin the 
data in x to extract a set of points for RdpixjQ 2 )). We combine the data points by 
taking the statistics- weighted average of the individual points, and the final systematic 
uncertainty is taken to be the statistics- weighted average of the individual systematic 
uncertainties. 

4. Nuclear Target Structure 

The current tensor F^ v of the target, 

F^(P,q) = Tv{T^(q) PP ) , T^{q) = I»{Q)I v {q) , (12) 
is a functional of the current-density operator 

/"(a) = e lP - x I»{0)e- lP - x , I v (q) = J d 4 qe^ 1 P \x) , (13) 

and the target density of a nucleus with spin j 

pp = 2fri E \°> p )( p >°\- ( 14 ) 

Lorentz transformations, A, are represented by unitary operators U(A). The 
current density operator and the target density transform according to 

U\K)T^(q)U{K)=k%,h v ^'"'(A" 1 ?) , 

U\N) PP U{K) = Pap . (15) 

The wave function representing the target state which generates that density may 
or may not be "manifestly" covariant. If the covariance is not "manifest" there 
always exist transformations to "manifestly" covariant representations, but there is no 
substantial significance to "manifest" covariance. Representations of target structure 
that emphasize the importance of "manifest" covariance intend to emphasize features 
of quantum field theory which affect the representation of the target structure. 

Assuming the only constituents of the deuteron are a proton and a neutron, the 
deuteron states \a, P) may be represented by functions ^p.M d .a(Pp, \,Pn, ^n), of 
the null- plane momenta and the null-plane spins. The null- plane spins are invariant 
under the kinematic Lorentz transformations that leave the null-plane invariant. These 
functions are eigenfunctions of an invariant mass operator, M and the four-momentum 
operator P: 

M 2 + P 2 

P±=P P ±+ Pn ±, P+=p++p+, P- = p+ ± . (16) 

The associated representation of the current tensor operator, 

(Pn> ^ ; Pp! ^'p\^ aa {Q)\Ppt \> Pm ^n), is the sum of the impulse term, 

(Pn> ^ri> Pp) ^pl-^"imp( < Z)|Ppi ^P) P"; ^n) 

= <5( P ;-p„)^ ;Atl (p;,A;i^( g )| Pj) ,A p ) 
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+ 5(p' p -p p )5 x , p , Xp {p' n ,\' n \F™(q)\p n ,\ n ) , (17) 

and the remainder, F^ig), which vanishes for large separation of the nucleons. The 
impulse term is essential for the relation of nucleon and deuteron structure functions. 
The representation oi F°^{q) is clearly model dependent. Any unitary transformation 
of p p which does not affect the N-N data [27] modifies Tf^t (l) an d thus anv conclusion 
relating nucleon and deuteron structure functions. For a = + and 2, the impulse 
assumption Tf^ = is consistent with the requirements of Poincare covariance and 
current conservation. 

It follows from the relation ([8]) and the impulse assumption that the deuteron 
structure function is related to the proton and neutron structure functions, 

F 2d {v,Q 2 ) = \(F 2p (y 1 Q 2 ) + F 2n {v 1 Q 2 )) , 

F 2N {v,Q 2 ) = f^r f d 2 P± ^p P (p)—F 2N (v N ,Q 2 ) , (18) 

J P + J P+ VN 

where vn = — (px ■Q)/m. No deep inelastic approximations are involved and there are 
no target mass effects. Target mass effects are artifacts of inconvenient representations 
of the deuteron states and the current density operators [28] . 

Nucleon momentum densities pp(p) obtained from a deuteron wave function 
satisfy the normalization condition 



/ 1£ / rfW(P) = i 



(19) 



With the definitions z = 2p + / P + and k± = pj_ — ^zP± the density p(z, k±) = pp(p) 
satisfies 

J'^ Jd 2 k ± p(z,k ± ) = l. (20) 

Conventional nuclear ground-state wave functions are eigenfunctions of a mass 
operator which are easily related to eigenfunctions of the four- momentum operator. 
In particular a deuteron wave function which satisfies 

(V + fc 2 /m) ip d =<p d E d , (21) 

where V is any "realistic" nucleon-nucleon potential, and is an eigenfunction of an 
invariant two-nucleon mass operator, 

(4cu 2 + 4mV)Mk) = Vd{k) M 2 d (l - E 2 /M 2 ) , (22) 



with to = \/k 2 + m 2 and 

Vda(k,ai,a 2 ) = |, <Ti, cr 2 |l, cr) 

+ Y / ^(\k\)Y 2 n (k)(l^,aua 2 \l,a')(2,l,n,a'\l,a) . (23) 
m,cr' 

The ratio E 2 d /Mj w 10" 6 is negligible. The invariant momentum density is 

p{k) = i- (| U (|fc|)| 2 + \w{\k\)\ 2 ) , J d 3 kp(k) = 1 . (24) 
The density p{z, fcjj = p(z, —k±) = p(2 — z, k±) is related to p(k) = p(—k) by 

r+oo 

p(z,kx)=z dk z p(k)5(z - 1 - k z /uo) . (25) 
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With these variable transformations the convolution relation of the structure 
functions, (fT5)) takes the form 

F 2N {v,Q 2 ) = f dz[d 2 k ± 1 —p( z ,k±) 

Jo J z- An 

x F 2N (v N , Q 2 ) 9 (mv N ~ \Q 2 ) , (26) 

with 

a k± Q± 2fc-Q ± 

A N = = — —z — x , v N = v(z - An) ■ (27) 

vnv Q z 

The Melosh rotations which relate null-plane spins to canonical spins do not 
affect these relations. Whether in practice the structure functions are parameterized 
as functions of v and Q 2 , x and Q 2 , or £ and Q 2 , is a matter of convenience. For 
interpolations and averaging of nucleon data, F 2 (£ i ,Q 2 ) has been used effectively. In 
the convolution, there is no particular justification for the use of of £, and a; is a 
more convenient choice. The kinematic constraint on the v and Q 2 dependence of the 
structure functions implies < x < 1 for the nucleons and < x < 2 for the deuteron, 
making it convenient to represent structure functions as functions F 2 (x,Q 2 ). In that 
representation the convolution relation, (f26|) . takes the following form: 

F 2N (x,Q 2 ) = [d 3 k—^ P (\k\)F 2N {x N ,Q 2 )0(l-x N ) , (28) 

J z(k)-A N 

where 

z(k) = I + k z /uj , x N = ^r— . (29) 

z - A N 

Note that in the case of the deuteron, k p — —k n , so z p _ n = 1 + fe^ ,n /u) = 1 ± k^/uj and 
A p = -A n . 

The main contributions are from values of z in the neighborhood of one (between 
0.8 and 1.2), and An is small compared to one, but not negligible. In the limit of 
An — > 0, the relation (f2"8"|) reduces to the familiar deep- inelastic convolution: 

F 2N (x,Q 2 )= J dzf{z)F 2N (- z ,Q 2 ) 6{z - x) , 

f{z) = [ d 2 k ±P (z,k ± ) . (30) 



5. Extraction of F 2n from Proton and Deuteron Data 



5.1. Extraction procedure 

Given the fit for F 2p from [14] . and the deuteron momentum density calculated from 
the CD-Bonn potential [29] , we use (f2"8")l to calculate F 2p . With a parameterization of 
F 2n we obtain F 2n and thus values of F 2 d, which are compared to the data. Because 
the ratio Rd P = F 2c i/F 2p is measured with smaller uncertainties than the absolute 
structure functions, we compare the calculated values of Rd p with the measured values 
to fit the parameterization of the neutron structure function. We parameterize R np at 
Q 2 = 12 GeV 2 and then use MINUIT [3D] to find the fit to R np that best reproduces 
the measured deuteron to proton ratio. 

At fixed Q 2 it is a matter of convenience whether one parameterizes R np — 
F 2n /F 2p as a function of x or £. We chose to parameterize R np as a function of £, with 
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the expectation that F2 n (£,Q 2 ) will have a smaller Q 2 dependence than F 2n {x,Q 2 ), 
as is the case for the proton and deuteron. We parameterize the R np (£,,Qo) using the 
following form, 

Rnp(0 = (Pi + P2O + P3 CXp (~p 4 f ) 

+p 5 exp(-p 6 (l -£)) +p 7 [max(0,C-p 8 )] 2 , (31) 

where the linear and first exponential terms are the dominant pieces at low x, and the 
quadratic and exponential terms provide flexibility at large x. At x = 0, R np f=a P1+P3, 
as the other terms are negligible. The initial fits gave values of R np consistent with 
unity at x = 0, and so we apply the constraint that the ratio should go to one, and 
take ps = 1 — p\ for all subsequent analysis. The last two terms yield deviation from 
the linear behavior at large x; by including two terms, we allow the possibility of 
significant cancellation between these terms, yielding flexibility for a modification of 
the shape at intermediate x values while still reproducing the data at large x. 

1.00 

0.95 

0.90 

r^f 0.85 

0.80 

0.75 

0.70 

0.65 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

X 

Figure 1. Extracted F^dlFlv f° r Q 2 =12 GeV 2 (points), along with the 
calculation taking R np from the best fit using the parameters in table [T] 



Table 1. Fit parameters for the best fit to R np using the parameterization from 
§31) . Note that p3 is not varied in the fit, and is taken to be (1 — pi). 



Parameter 


value 


Parameter 


value 


Pi 


0.816 


P5 


-0.034 


P2 


-0.661 


P6 


8.714 


(ps) 


0.184 


P7 


-0.072 


P4 


5.509 


PS 


0.450 



Figure [T] shows the Rd P data interpolated to Q\=Yl GeV 2 . The uncertainties 
shown are the combined statistical and uncorrelated systematic uncertainties. There 
is also an estimated 1% normalization uncertainty [16j . The curve in figure [1] shows 
the calculated Rd P (x,Qg) for the best fit to R np (the parameters of the fit are given 
in table [T]) . Because the uncertainties are almost evenly split between statistical and 
correlated systematic except at the largest x values, we expect a reduced chi-squared 
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Figure 2. Extracted Fi n jFi v for Q 2 =12 GeV 2 . The data points represent values 
extracted from Rd p measurements according to H33jl while the curve is given by 
(13 1 I t with the parameters in table [T] 



of less then unity. Including the contribution from the normalization paremeters, the 
fit gives x 2 = 14.6 for 17 degrees of freedom (27 data points, 7 parameters in the 
^-dependence, and 3 normalization factors). 

Figure[2]shows the fit for R np , along with the values of R np extracted from the Rdp 
measurements. In terms of the smearing ratios, Sn(x, Q 2 ) — F2n{x, Q 2 )/ F2n(x, Q 2 ), 
we can write (TT51) as: 

Fid — \{S P F2 V + S n F2n) , (32) 

which allows us to extract R np using the simple expression: 

R np = (2-Rdp — S p )/S n - (33) 

The advantage of extracting individual R np points using these smearing ratios, 
calculated from the extracted F2 n (x, Qq), is that it provides an estimate of the 
uncertainty in R np coming from the uncertainty in the Rd p measurements. While this 
is convenient, and consistent with the way previous extractions have been presented, it 
does not take into account that a given value of Rd p (x, Q 2 ) depends on the proton and 
neutron structure functions over a range in x. This will be discussed in section [5.21 



5.2. Systematic uncertainties 

The extracted R np (x) values shown in figure [2] include the uncertainties associated 
with the Rd p data, but not the uncertainties associated with the deuteron model and 
the impulse assumption or systematic uncertainties associated with the extraction 
procedure. These uncertainties tend to be important at large x, and will tend to 
have a highly correlated impact on the high- a; results, and thus the extrapolation 
to x = 1. Because the systematic uncertainties were observed to grow with x, and 
yield highly correlated corrections among the high- a: points, we make estimates of the 
uncertainties and then fit these to an exponential, to provide a simple parameterization 
of the correlated systematic error band. 
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The largest systematic at all x values was due to the 1% overall normalization 
uncertainty on the global Rd p measurements. Shifting the deuterium ratios by 1% 
yields an overall absolute shift in R np of approximately 0.014, plus an additional 
contribution at large x values. The large x contribution is well described by 5R np = 
0.079 exp (-12(1 -at)). 

The result is also sensitive to the parameterization of the proton structure function 
at large x. We vary the high- a; behavior of F2 P (x, Q 2 ), introducing a 10% change in the 
falloff at large x, which is large enough that the F2 P fit starts to become inconsistent 
with the SLAC measurements. This yields a change of 5R np = 0.06 exp (—12(1 — x)), 
where we have used the same x dependence as for the normalization systematic, as 
the behavior is very similar. 

Q t Q r — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — 




o.o r 1 

0.4 0.5 0.6 0.7 0.8 0.9 1.0 

X 

Figure 3. Extracted value of R n p, along with the best fit (solid line) and 
systematic uncertainty bands (dotted lines). 

Finally, the uncertainties shown on the R np data point come only from the 
statistical and systematic uncertainties on the Rd P ] there is no contribution associated 
with the uncertainty in the factor S n , which depends on the extracted F^n- To account 
for this, we vary the parameters in the fit to R np and calculate Rd P - We compare the 
extracted Rd P to the data, and take the range of fits that increase the total \ 2 by one 
as the range of good R np fits. For each of these, we recalculate the smearing ratio, 
S n , and determine the impact on the extracted R np data points. Again, this effect is 
largest at large x, and can be parameterized as 5R np = 0.05 exp (—12(1 — x)). 

Combining these systematic uncertainties in quadrature, we find a overall 
uncertainty of 0.014, and a correlated uncertainty that grows with x of the form 
SR np — 0.11 exp (—12(1 — x)). Figure [3] shows the extracted values of R np at 
Q 2 = 12 GeV 2 , along with the error band associated with the systematic uncertainties. 
Our extracted values, rebinned by a factor of two and including both statistical and 
systematic uncertainties are given in Tabled 
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Table 2. Extracted values and uncertainties for R dp (x) and R n p(x) at 
Q 2 =12 GeV 2 . Errors labeled "stat" are the combined statistical and uncorrelated 
systematic uncertainties. Uncertainties labeled "sys" are the correlated systematic 
uncertainties coming from the relative normalization between different data sets 
and are correlated. The final column is the uncertainty in the extracion of R np due 
to the systematic uncertainties in the extraction procedure discussed in Sec. 15.21 
Not included is the overall 1% normalization uncertainty in R dp which corresponds 
to an overall shift of ±0.014 in R np . 



X 


Rdp 


5R dp at 




Rnp 


0H np 


SRnp° X 


SR-np 1 


0.030 


0.976 


0.006 


0.005 


0.954 


0.013 


0.010 


0.000 


0.050 


0.964 


0.006 


0.005 


0.929 


0.013 


0.010 


0.000 


0.079 


0.938 


0.006 


0.005 


0.877 


0.012 


0.011 


0.000 


0.105 


0.927 


0.008 


0.008 


0.855 


0.016 


0.015 


0.000 


0.140 


0.906 


0.008 


0.007 


0.815 


0.015 


0.014 


0.000 


0.180 


0.875 


0.008 


0.007 


0.753 


0.015 


0.014 


0.000 


0.226 


0.855 


0.008 


0.007 


0.716 


0.015 


0.014 


0.000 


0.275 


0.831 


0.008 


0.007 


0.670 


0.015 


0.014 


0.000 


0.324 


0.819 


0.010 


0.006 


0.646 


0.019 


0.011 


0.000 


0.353 


0.801 


0.007 


0.007 


0.610 


0.014 


0.013 


0.000 


0.425 


0.775 


0.006 


0.005 


0.559 


0.012 


0.009 


0.000 


0.458 


0.766 


0.008 


0.007 


0.539 


0.016 


0.014 


0.000 


0.519 


0.747 


0.006 


0.005 


0.497 


0.012 


0.010 


0.000 


0.563 


0.725 


0.007 


0.006 


0.447 


0.015 


0.013 


0.001 


0.600 


0.718 


0.006 


0.005 


0.428 


0.012 


0.010 


0.001 


0.625 


0.712 


0.008 


0.006 


0.410 


0.016 


0.012 


0.001 


0.660 


0.709 


0.007 


0.006 


0.391 


0.014 


0.012 


0.002 


0.690 


0.713 


0.006 


0.005 


0.383 


0.012 


0.009 


0.003 


0.727 


0.712 


0.008 


0.006 


0.357 


0.015 


0.012 


0.004 


0.756 


0.707 


0.007 


0.004 


0.323 


0.012 


0.008 


0.006 


0.780 


0.716 


0.009 


0.006 


0.311 


0.015 


0.011 


0.008 


0.805 


0.734 


0.009 


0.006 


0.304 


0.016 


0.010 


0.011 


0.826 


0.736 


0.011 


0.005 


0.267 


0.018 


0.009 


0.014 


0.861 


0.815 


0.025 


0.008 


0.297 


0.037 


0.012 


0.021 


0.870 


0.841 


0.020 


0.007 


0.304 


0.028 


0.011 


0.023 


0.887 


0.821 


0.027 


0.009 


0.212 


0.036 


0.012 


0.028 


0.901 


0.833 


0.038 


0.009 


0.169 


0.047 


0.011 


0.034 


tAdditional 1% scale uncertainty is not shown. 


f Additional scale uncertainty of 0.014 is not shown. 



5.3. Q 2 dependence 

While the analysis is performed at fixed Q 2 , we fit R np as a function of £, with the 
expectation that the Q 2 dependence will be minimal when taken at fixed £, as is 
observed for the proton and deuteron structure functions. By assuming that this 
parameterization is a reasonable representation of the Q 2 dependence of the ratio, we 
can use this and the Q 2 -dependent fit to the proton structure function |14j to examine 
the Q 2 dependence of the calculated Rd p . 

Figure |4] shows Rd P calculated at Q 2 =6 (bottom curve), 12, and 20 GeV 2 (top 
curve). There is a significant Q 2 dependence at large x, coming mainly from the Q 2 
dependence of F-i-p. We compare this with the data, now taken in smaller Q 2 ranges (4- 
8, 8-16, and 16-32 GeV 2 ), interpolated to the same Q 2 values as the curves. While 
the data at large x is limited in the Q 2 range, there is a clear Q 2 dependence to 
the measured ratios which is at least in qualitative agreement with the result of the 
calculation based on the extracted R np (0- 
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Figure 4. Ratio of deuteron to proton structure function for Q 2 values from 
4-8 GeV 2 (black circles), 8-16 GeV 2 (red crosses), and 16 - 32 GeV 2 (blue 
diamonds), compared with calculations at Q 2 =6,12, and 20 GeV 2 , using the 7- 
parameter fit to R n p{^) shown in figure [2] 



6. Discussion of the Results 



In the deep-inelastic limit, x 2 -C Q 2 /4m 2 , structure functions are related to parton 
distributions d(x) and u{x). In that limit, 

d(x,Q 2 ) m np (x,Q 2 )-l 

u{x,Q 2 ) ~ A-R np {x,Q 2 ) ' 1 ' 

plus small contributions from heavier quarks. For Q 2 =12 GeV 2 this is expected to 
be valid for x <, 0.5. Our extraction of R np (x) is consistent with that calculated from 
the CTEQ6L [31 parton distributions at all x values, although both extractions have 
large uncertainties at the larger x values, where the simple connection to the quark 
distributions is expected to break down. 

At large x, the quantitative features of the extracted neutron structure function 
are limited by the availability of high precision deuteron and proton data. Existing 
measurements provide a reasonable basis for the extraction of R np with high precision 
up to x = 0.85, but no significant constraint for x > 0.9. It follows from the proton 
and deuteron data, and the deuteron structure assumptions we specified, that the 
neutron to proton ratio is small for large x, although there is little to constraint the 
extrapolation to x = 1. 

In figure[5]we compare the R np ratio obtained in section[5]to previously published 
values of Fi n jFi p . The solid squares are the results of Whitlow et al. [16], who use 
the DIS convolution and the deuteron wavefunction of the Paris potential |33j . The 
hollow squares are the results of [32], which emphasizes the off-shell effects of the 
Buck-Gross [34] spectator representation of the deuteron. 

There are a few difficulties in comparing these previous results to our extraction. 
As mentioned in section [5721 the previous extractions compared in figure [5] neglect the 
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Figure 5. Extracted Fin/Pip for Q 2 =12 GeV 2 (solid circles) compared to 
previous extractions. The hollow squares are points from |32j . while the other 
points are from I16| , using a DIS convolution (solid squares) and using a density- 
dependent extrapolation of the EMC effect (hollow circles). 




Figure 6. The smearing ratio, Sp, calculated at several Q 2 values. The solid 
line is 12 GeV 2 , while the dotted lines correspond to the Q 2 values of the R^ p 
points extracted by Whitlow [16| . and used in previous analysis. The lowest curve 
corresponds to Q 2 = 4.7 GeV 2 , while the upper curve is 23.6 GeV 2 . The points 
are placed on the curve associated with the average Q 2 for that x value, which is 
significantly different from the 12 GeV 2 value for the large x points. 



fact that the uncertainty in the fits to F 2p and the extracted F% n affect the smearing 
ratios S p and S n (|32|) , and thus their total uncertainties are close to our statistical 
uncertainties, as most of the additional systematics are neglected. In addition, the x 
binning is very coarse. For x > 0.7, there is a significant ir-dependence to Rd P - If the 
average x of the data in one of the large x bins is shifted from the center of the bin, 
placing the point at the central x value will affect the results. 
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More important is the strong dependence of the smearing ratios on Q 2 , combined 
with the strong variation of the average Q 2 value used for each x bin which affects the 
results compared in figure These extractions used Rd p values from [16] , where each 
x value has a different average Q 2 . The low x points have (Q 2 ) values below 8 GeV 2 , 
while the highest x point has (Q 2 )=23.6 GeV 2 . The extraction of [32] calculates S p 
and S n at Q 2 =12 GeV 2 , and uses this for all x values. However, the smearing functions 
have a strong dependence on Q 2 at large x values, as shown for S p in Figure O For 
x = 0.85, the difference between S p at 23.6 GeV 2 at 12 GeV 2 (the solid line) is 12%, 
and so the proton contribution to Rd P should be approximately 0.12 higher, making 
R np roughly 0.12 lower. A precise comparison of the results obtained using these 
different models for the deuteron structure contributions will require a self-consistent 
treatment of the data in the extraction of R np - 

If one assumes [35] that the ratio 2F2d/{F%p + ^2n) can be obtained by scaling 
the EMC ratio, F 2 A/F2d, then R np can be extracted from F 2 dlFi P data without 
calculating smearing ratios S p and S n . Such an extraction was performed in [16j . 
yielding the largest values (hollow circles) for the ratio in figure [5] This procedure 
involves the questionable assumption that the details of the nucleon momentum 
distribution in the deuteron are not important j36l I37j . Further, it entirely neglects 
the important difference between S p and S n by extracting a universal smearing ratio 
from the extrapolation of the EMC ratio. 

No clear conclusions can be drawn from comparisons of the results of different 
extractions unless they take consistent proton and deuteron data as input and properly 
treat the Q 2 dependence of the data. Given proton and deuteron data at fixed Q 2 , the 
effect of different deuteron structure assumptions should be investigated. It will be 
important not to confuse different structure assumptions with different representations 
of the same deuteron structure. Here we assumed that the deuteron state is represented 
by a vector in the two-nucleon Hilbert space (tensor product of two single-nucleon 
Hilbert spaces), together with the impulse current tensor. The deuteron wave function 
which implements these assumptions is an eigenfunction of a mass operator which 
implements one-pion exchange and two-pion exchange with some parameterization of 
the short range features which affects the high momentum features of the momentum 
density p(k). This is implemented differently by different potentials which produce 
high precision fits to nuclcon-nucleon scattering data. 

Relations to larger nuclei are important. The impulse assumption with null- plane 
kinematics applied here to the deuteron is easily applied to to larger nuclei, where 
eigenfunctions of realistic mass operators exist. Ratios of 4 He/ 2 H will test the the 
smearing effects on the average nucleon structure functions while 3 He/ 2 H ratios will 
provide an additional test of both the nuclear effects and the extracted the neutron 
structure function New measurements of these ratios at large x will be available soon 
from a recent Jefferson Lab measurement [38] . 

7. Conclusion 

The relativistic quantum theory of nuclear structure provides the basis for reliable 
extraction of neutron structure functions from proton and deuteron structure 
functions. We extract the neutron structure function using an impulse approximation 
and an effective two-nucleon mass operator. No deep inelastic approximations or quark 
structure assumptions are involved. We find that sufficiently accurate interpolations of 
the proton and deuteron data to fixed Q 2 are essential, and quantitative comparisons 
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to other extraction procedures are not possible until these extractions include a proper 
treatment of the Q 2 dependence of the data. 

We obtain a precise extraction of the neutron structure function up to x w 0.85, 
within the context of our dcuteron structure approximation. For x > 0.9, there are 
no data to constrain the calculation, and any extrapolation of the results to larger x 
is unreliable. More data, in particular at large x values, will be essential to extending 
the extraction. 
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